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The method of iterated resolvents is used to obtain an effective 
Hamihonian for neighbouring qubits in the Kane sohd state quantum 
computer. In contrast to the adiabatic gate processes inherent in the 
Kane proposal we show that free evolution of the qubit-qubit system, 
as generated by this effective Hamiltonian, combined with single qubit 
operations, is sufficient to produce a controlled-NOT (c-NOT) gate. 
Thus the usual set of universal gates can be obtained on the Kane 
quantum computer without the need for adiabatic switching of the 
controllable parameters as prescribed by Kane Both the fidelity 
and gate time of this non-adiabatic c-NOT gate are determined by 
numerical simulation. 
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I. INTRODUCTION 



In the Kane proposal for a solid state quantum computer a series of spin | ^^P 
nuclei in a silicon substrate are used as qubits The interaction between these 
qubits is mediated by valence electrons weakly bound to the nuclei, such that at 
energies much lower than these electrons binding energy the system Hamiltonian is 
given by 

+ AiCTie-Cfin + ^20^2e-'?2n + Jaie.(?2e- (1) 

The interaction strengths Ai , A2 as well as J are controllable by means of voltage 
biases applied to "A" and "J" gates respectively [00. In the case where J = 0, 
the energy splitting between the electron spin states is approximately 1600 times 
larger than the splitting of the nuclear states. Thus at temperatures low compared 
to these energy splittings the nuclear spin states, in the electron ground state, can 
be manipulated without significantly altering the electron states. This gives a basis 
in which quantum computation can be performed, the reduced Hilbert space of the 
electron ground state. The quantum computing basis is then: 

UiOO),|ii01),|iilO),|iill). (2) 

Here | |) denotes an electron in the spin down state and |0), |1) denote a nucleus in 
the spin up and down state respectively. The electron spin down state remains the 
ground state as long as J < although the energy difference between this state 
and the odd superposition of electrons reduces as J increases, thus in operating our 
quantum computer we want to keep J below this limit. 

A requirement of a universal quantum computer is that it be able to implement 
a universal set of gates. One such universal set comprises of the set of single qubit 
rotations and the controlled-NOT (c-NOT) gate The currently proposed im- 

plementation of the c-NOT in the Kane quantum computer relies on an adiabatic 
switching on of the A and J couplings to produce a unique energy splitting between 
states which can then be swapped using a Rabi type flipping induced by a transverse 
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magnetic field rotating at a frequency sucli tliat it is in resonance witli tlie desired 
transition . Tliis adiabatic switcfiing process lias been studied in detail [|TT| and 
can produce a c-NOT gate with an error of e ~ 10^^ in a time tc-NOT = 26/is. In 
this article we introduce a c-NOT gate that does not rely on adiabatic switching. 
We derive an effective Hamiltonian that describes the interaction of the qubits in 
the subspace of the electron ground state using the method of iterated resolvents. 
This effective Hamiltonian is then used to derive a c-NOT gate that relies only on 
the evolution of the qubits that it generates, in combination with the single qubit 
operations. We then use this fact to argue that the family of gates consisting of the 
free evolution of neighbouring qubits for pre determined times, is a suitable alter- 
native to the adiabatically constructed c-NOT gate for use in the implementation 
of quantum algorithms. 



II. THE EFFECTIVE HAMILTONIAN 



To understand the dynamics of the qubits in the computational subspace, it is 
useful to calculate an effective Hamiltonian that describes the action of Eq (|I|) in 
this reduced bases. To do this we use the method of iterated resolvents 0, which 
has been used successfully to calculate effective Hamiltonians in reduced bases in 
QCD. The method of iterated resolvents involves the systematic reduction of the 
dimension of the system to that of the subspace of interest. In this case a reduction 
of the full 16 X 16 system Hamiltonian, i^ig, to a 4 x 4 effective Hamiltonian, H4. 
We begin by writing the eigenvalue problem for the complete Hamiltonian 

16 

Y.mM{m = Em. (3) 

Let us divide the rows and columns of Hi^ into separate subspaces, take the first 15 
rows to be the first subspace, call it P, and the remaining subspace we call Q. Thus 
Eq(|) can be expressed as 

f2{t\H^,\j){m + {^\H^e\16){lQm = E{z\^il), (4) 
i=i 

which can be written in the block matrix form 
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(p|iji6|p)(p|^) + (p|i/i6|g)(Q|^) = (5) 

{Q\H,,\P){P\^) + {Q\HM{Q\^) = E{Q\^). (6) 

Because the eigenvalue E is in general unknown, it is replaced with a free parameter 
uj to be determined later. Thus if the matrix {Q\uj — Hiq\Q) can be inverted, the Q 
space wave function can be expressed in terms of the P space wave function. 

{Ql^iuj)) = Gq{uj){Q\H,,\P){P\^), (7) 

where we have written the resolvent 

Substituting Eq(^ into Eq(|^) gives 

{{P\H,e\P) + {P\HMGQ{u){Q\H,e\P)){Pm = PM(P|^), (9) 
which defines an eigenvalue equation in the P space 

{P\Hi,{u;)\P){P\^iu)) = Eiu){P\^{u)), (10) 
which in turn defines for an effective 15 Hamiltonian in the P space 

ifisH = ^16 + H,e\Q)GQ{u){Q\H^e. (11) 

The method of iterated resolvents calls for this procedure to be repeated until an 
effective 4x4 Hamiltonian for the computational sub-space is produced. Finally it is 
necessary to solve the fixed point equation (4|if4(w)|4) = w. This equation produces 
many solutions, the correct one to choose is the one that yields an eigenspectrum 
for Hi that is as close as possible to the four lowest eigenvalues of the complete 
Hamiltonian. This procedure was completed numerically for parameters similar to 
those used in the adiabatic gate Ai = A2 = 1.683 and J = 600 in units of gnfJ-nB. 
This yielded an effective Hamiltonian of the form 

H,s = A((tJ^4 + af4) + Qalal + A(cr^ + a'^) + F, (12) 

where it is understood that the subscripts denote the first and second nucleus. 
The values of the parameters were A = 2.3723 x 10"^ 6 = -1.4645 x 10"^A = 
2.6871, F = 5.3578. 
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III. CONSTRUCTING A C-NOT GATE 



The generator of the c-NOT operator, f/c_NOT = exp[— iC], is given by 

G=^-{l-al-al + ala^,). (13) 

All the terms on the right hand side commute, and so the c-NOT operator can be 
written as 

[/c-NOT = exp[-^]exp[^ai"]exp[^(T2^exp[-^aiV^]. (14) 

Here, reading from left to right, the first factor just represents a phase factor, this 
is unimportant as all that is required is a c-NOT operation up to an overall phase. 
The second factor is a rotation, this can be realized by a combination of the 
single particle rotations common to NMR theory [§,|l3 . 

expi^a^ = exp[^(T"]exp[-^(T^]exp[-^cr"]. (15) 

The third factor is just another standard single particle operation. The fourth factor 
requires a combination of pulses to give the required evolution. Given an effective 
Hamiltonian of the form Eq ([l^) we can use standard refocusing techniques standard 
to construct this evolution of the qubit-qubit system. The c-number factor in the 
effective Hamiltonian commutes with all other terms and simply leads to an overall 
phase in the evolution, it can thus be ignored in the calculations and included at 
the end. Our first step then is to refocus out the Zeeman evolution of the nuclei, 
using rf pulses targeted at both nuclei simultaneously: 

ZTT — xty 

exp[-(a^ + (^2)] ^M^iHeff - r)] X 
exp[-y (err + a^,)] exp[— (ife// - T)] 

= exp[^(A(ara2^' + a^f ) + Qala^)]. (16) 

The next step is to refocus out the erf erf and parts of the evolution: 

777- — 
exp[-at] exp[— (A«a2^- + afaf) + Oa^V^^)] x 

?7r — 7f 

exp[-ya[] exp[— (A(ara| + afaf) + Qa^a^,)] 

= expl^Aa^a^l. (17) 



Again targeting the first spin only, we can obtain 

%1X — it — ZTT — it 

exp[-ar] exp[— Aa^a2"] exp[^a^] = exp[— AaiV2"]; (18) 

this evolution for a time t = 7r/i/(4A) gives the required operator in Eq (|1^). Sub- 
stituting these pulses into Eq(|l4D we find that the c-NOT operation can be written 
as, 

Uc_NOT = exp[- — (1 - — )]exp[— cri]exp[-— cjf]exp[-— 
X exp[ya2lexp[^af]exp[ya[]exp[y (or^= + a^)] 

X F\ — T- exp (o-f + (Tn) F — - exp — (To F — - 

X F[^]exp[-^ar], (19) 

where F[t\ = exp[iTTtHeff /h] denotes free evolution generated by -f^e// over a time t. 
The first term is simply a phase correction, it cannot be implemented physically but 
tells us by what overall phase the composite operator must be corrected to produce 
exact c-NOT evolution. 



IV. GATE TIME AND FIDELITY 

It is first necessary to check that Eq(|l^) does give an accurate description of 
the two qubit effective Hamiltonian. This was done by numerically solving the 
Schrodinger equation for the time development of the entire electron-nuclear system 
using the Hamiltonian Eq(||) over a time of several microseconds, and comparing 
the evolution of the qubits to that predicted by the effective Hamiltonian over the 
same period. It was found that final states agreed with an error probability of 10^^. 

Let us now calculate the time this non-adiabatic gate takes to execute, using 
similar values for the operating parameters to those prescribed for the adiabatic 
c-NOT gate. The construction requires that the system evolves freely for a total 
time t = 7r/(4A) ^ 3fis. In addition to this free evolution, we also have the single 
particle rotations, the time scale of which is set by the time it takes to implement 
a TT rotation, r ^ 22/is. The combination of single particle operations necessary for 
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this implementation takes a total of t ~ 77/is to execute, taking the total gate time 
to tc-NOT ~ 80/is, compared with 26/is for the adiabatic implementation. 

This gate was simulated numerically and it was found that the operation is indeed 
that of a c-NOT gate, with an error probability of e = 4 x 10~^. This is slightly 
outside current estimates of the error tolerance of a QC using error correcting codes, 
which vary between e < 10^^ — 10^^ [P,p|,|T2|. These error tolerances, most authors 
admit, are probably pessimistic, and are based on some very general assumptions 
about the type of error process and the architecture of the computer. It is possible 
that by tailoring an error correcting code to a specific problem, the error tolerance 
may be more forgiving and the error probability of this gate may fit well within the 
new bound. 

The problem still remains however, that this non-adiabatic c-NOT gate is both 
slower, and of lower fidelity than the adiabatic implementation. Both these facts 
can in part, be attributed to the large number of single qubit operations performed 
in the c-NOT implementation. Note that the period of "free evolution" required 
to implement the c-NOT gate is only about 3/is. During this free evolution is the 
only time that there can possibly be information flow between the two qubits, in 
contrast to the adiabatic case in which information flow occurs over the entire gate 
time of 26/is. With this in mind, we can consider free evolution for this time to be 
an elementary two-qubit gate that is capable of transferring the same quantity of 
quantum information between two qubits as is a c-NOT gate, in just over 10% of the 
time. It therefore seems reasonable to speculate that it would be possible to recast 
any large scale quantum algorithm, into a series of single qubit rotations, and two- 
qubit "free evolutions" without the need to explicitly construct a c-NOT gate, and 
that this recipe for constructing the algorithm may be faster than one which relies 
on an explicit adiabatic implementation of the c-NOT gate, for the Kane quantum 
computer. 
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